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PREFACE 


This  report  describes  a  method  for  estimating  the  flov  rate  through  the  gaps 
of  offshore  segmented  breakwaters  caused  by  wave  overtopping.  Factors  that  can 
be  investigated  using  this  method  are  the  influence  of  breakwater  freeboard, 
wave  height  and  period,  breakwater  length  and  spacing,  number  of  breakwaters, 
distance  offshore,  water  depth  at  the  breakwater,  and  shore  attachment  on  the 
flow  rate.  Other  wave  effects  on  hydraulics,  such  as  diffraction,  refraction, 
reflection,  and  wave-current  interactions,  have  not  been  considered.  The  work 
was  carried  out  under  the  offshore  breakwaters  for  shore  stabilization  and 
evaluation  of  shore  protection  structures  programs  of  the  U.S.  Army  Coastal 
Engineering  Research  Center  (CERC). 

The  report  was  prepared  by  William  N.  Seelig  and  Dr.  Tcdd  L.  Walton,  Jr., 
Hydraulic  Engineers,  under  the  general  supervision  of  Dr.  R,M.  Sorensen,  Chief, 
Coastal  Processes  and  Structures  Branch  and  Dr.  J.R.  Weggel,  Chief,  Evaluation 
Branch . 

Comments  on  this  publication  are  invited. 


Approved  for  publication  in  accordance  with  Public  Law  166,  79th  Congress, 
approved  31  July  1945,  as  supplemented  by  Public  Law  172,  88th  Congress, 
approved  7  November  1963. 


Colonel,  Corps  of  Engineers 
Commander  and  Director 
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CONVERSION  FACTORS,  U.S.  CUSTOMARY  TO  METRIC  (SI)  UNITS  OF  MEASUREMENT 


U.S.  customary  units 
metric  (SI)  units  as 

of  measurement 
follows: 

used 

in  this  report  can  be  converted  to 

Multiply 

by 

To  obtain 

inches 

25.4 

millimeters 

2.54 

centimeters 

square  inches 

6.452 

square  centimeters 

cubic  inches 

16.39 

cubic  centimeters 

feet 

30.48 

centimeters 

0.3048 

meters 

square  feet 

0.0929 

square  meters 

cubic  feet 

0.0283 

cubic  meters 

yards 

0.9144 

meters 

square  yards 

0.836 

square  meters 

cubic  yards 

0.7646 

cubic  meters 

miles 

1.6093 

kilometers 

square  miles 

259.0 

hectares 

knots 

1.852 

kilometers  per  hour 

acres 

0.4047 

hectares 

foot-pounds 

1.3558 

newton  meters 

millibars 

1.0197  x 

10"3 

kilograms  per  square  centimeter 

ounces 

28.35 

grams 

pounds 

453.6 

grams 

0.4536 

kilograms 

ton,  long 

1.0160 

metric  tons 

ton,  short 

0.9072 

metric  tons 

degrees  (angel) 

0.01745 

radians 

Fahrenheit  degrees 

5/9 

Celsius  degrees  or  Kelvins1 

xTo  obtain  Celsius  (C)  temperature  readings  from  Fahrenheit  (F)  readings, 
use  formula:  C  -  (5/9)  (F  -32). 

To  obtain  Kelvin  (K)  readings,  use  formula:  K  -  (5/9)  (F  -32)  +  273.15. 
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SYMBOLS  AMD  DEFINITIONS 


a,  b 


B 


d 

8 


F 

g 

H 


H' 

o 

H 

8 

h 


K 

L 

o 

*, 


in 

N 

P 

Q 

% 

«n 

% 

R 


cross-sec Clonal  flow  area  of  gap  between  breakwaters 
cross-sectional  flow  area  between  breakwaters  and  shoreline 
empirical  runup  coefficients 
breakwater  gap  width 

discharge  coefficient  for  breakwater  gaps 

discharge  coefficient  for  space  between  breakwaters  and  shoreline 
water  depth  at  the  toe  of  the  structure 
breakwater  freeboard  -  h  -  d 

s 

acceleration  due  to  gravity 
wave  height  at  the  structure 
equivalent  deepwater  wave  height 
significant  wave  height  at  the  structure 
structure  height 

difference  in  mean  water  levels  inside  and  outside  of  breakwaters 

a  dimensionless  parameter 

deepwater  wavelength 

breakwater  length 

log  e 

number  of  breakwaters 
ponding  level 

mean  net  flow  rate  through  a  breakwater  gap  or  Inlet 
net  inflow  by  overtopping  of  a  breakwater 
empirical  overtopping  parameter 

net  inflow  by  overtopping  per  unit  length  of  a  breakwater 
overtopping  per  unit  length  of  breakwater  with  no  return  flow 
runup 
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SYMBOLS  AND  DEFINITIONS — Continued 

vave  period 

mean  net  water  velocity  for  flow  through  a  gap  -  Q/Ac 
dimensionless  velocity 
empirical  overtopping  parameter 
breakwater  seaward-face  slope  angle 


surf  similarity  parameter 


ESTIMATION  OF  FLOW  THROUGH  OFFSHORE  BREAKWATER 
GAPS  GENERATED  BY  WAVE  OVERTOPPING 

by 

William  N .  Seelig  and  Todd  L.  Walton,  Jr. 

I.  INTRODUCTION 

Offshore  breakwaters  are  often  constructed  to  protect  harbors  or  eroding 
coasts  from  wave  action.  Wave  overtopping  of  segmented  breakwaters  generates 
a  seaward  flow  through  breakwater  gaps.  Low  discharges  and  low  water  velocities 
are  usually  desirable  in  gaps  between  breakwaters.  High  velocities  may  be  a 
hazard  to  navigation  or  swimmers,  and  large  net  exit  flows  through  the  gaps 
may  transport  sediment  out  of  the  breakwater  area.  This  report  Illustrates  a 
technique  for  predicting  the  net  discharges  and  mean  velocity  through  breakwater 
gaps  caused  by  overtopping. 

II.  HYDRAULICS  OF  DETACHED  BREAKWATERS 

Since  the  cost  of  a  breakwater  increases  with  the  height  of  the  structure, 
it  may  be  necessary  to  build  the  structure  to  allow  some  wave  overtopping  (Fig. 
1) .  A  moderate  overtopping  rate,  which  helps  maintain  water  quality  by  encour¬ 
aging  circulation,  may  also  be  desirable.  If  the  breakwater  is  impermeable  the 
volumetric  rate  of  overtopping  per  unit  length  by  breakwater,  qQ  (assuming 
no  return  flow  over  the  breakwater),  may  be  predicted  (Weggel,  1976;  U.S.  Army, 
Corps  of  Engineers,  Coastal  Engineering  Research  Center,  1977)  by 


0.1085 

PROTECTED 


9 


where 


g  *  the  acceleration  due  to  gravity 

Hq  =  the  equivalent  unrefracted  deepwater  wave  height 

R  -  the  vertical  height  of  runup  on  the  structure  if  the  breakwater 
were  high  enough  so  that  no  overtopping  occurred 

F  =  breakwater  freeboard  =  h  -  dg 

h  =  the  structure  crest  elevation 

dg  =  the  water  depth  at  the  toe  of  the  structure. 

Q*  and  a  are  empirical  overtopping  coefficients  found  in  the  Shore  Protection 
Manual  (SPM)  (U.S.  Army,  Corps  of  Engineers,  Coastal  Engineering  Research  Center, 
1977,  Ch.  7).  Sample  values  of  the  empirical  coefficients  for  an  impermeable 
riprap  structure  with  a  1  on  1.5  seaward  slope  are  shown  in  Figure  2. 

A  first  approximation  of  the  wave  runup  on  rubble-mound  breakwaters  may  be 
estimated  using  the  equation  of  Ahrens  and  McCartney  (1975): 


R  = 


(i  +  bO’  * 


tan  0 


(2) 


where  L0  is  deepwater  wavelength,  and  8  is  the  slope  angle  of  the  seaward 
face  of  the  breakwater;  the  empirical  coefficients  a  ■=  0.692  and  b  =  0.504  are 
recommended.  Note. — Other  methods  for  estimating  runup  on  various  structures 
may  be  found  in  the  SPM  (U.S.  Army,  Corps  of  Engineers,  Coastal  Engineering 
Research  Center,  1977)  and  Stoa  (1979). 

Equations  (1)  and  (2)  were  developed  using  monochromatic  wave  tests,  so 
they  should  be  used  for  swell  wave  conditions  where  the  wave  height  and  period 
from  one  wave  to  the  next  is  approximately  constant .  The  method  of  Ahrens 
(1977)  for  modifying  equations  (1)  and  (2)  for  irregular  waves  is  recommended 
for  irregular  waves  generated  by  nearby  storms.  The  irregular  wave  overtopping 
prediction  procedure  is  rather  complicated,  so  the  computer  program  BWFL0W2 
(CERC  program  number  752X6R1ANC)  is  recommended  for  irregular  waves.  This  pro¬ 
gram  is  available  '  i  the  Corps  of  Engineers  Computer  Library  at  U.S.  Army 
Waterways  Experiment  Station,  Vicksburg,  Mississippi.  Note  that  the  irregular 
wave  overtopping  method  tends  to  be  conservative  because  a  Rayleigh  wave  height 
distribution  is  assumed,  while  the  actual  distribution  may  be  truncated  due  to 
depth  or  steepness  limited  breaking. 

1.  Enclosed  Breakwater  Systems. 

If  the  breakwater  system  is  enclosed  on  either  end  by  impermeable  groins 
and  the  breakwater  has  no  gaps,  water  overtopping  the  breakwater  would  cause 
the  water  level  landward  of  the  breakwater,  h^,  to  rise.  Eventually,  the 
zone  landward  of  the  breakwater  would  fill  up  to  a  ponding  level  where  the  sea¬ 
ward  flow  of  water  over  the  breakwater  would  equal  inflow  and  the  net  flow,  qn, 
would  be  zero.  Diskin,  Vajda,  and  Amir  (1970)  tested  a  number  of  breakwaters 
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(ft  rt«'| 


at  various  water  depths,  wave  heights,  and  wave  periods  and  found  that  the 
ponding  level  (Fig.  3)  may  be  estimated  from 


|r  -  0.6*  *  5j) 


0.1  for  >  2.05 
Ho 


(3) 


where  F  is  the  breakwater  freeboard  defined  as  (h  -  ds).  The  maximum  ponding 
level  occurs  at  F/Hq  =0.7  because  this  is  the  point  just  before  the  seaward 
flow  over  the  top  of  the  breakwater  occurs.  Note  that  there  is  scatter  in  the 
data  used  to  develop  equation  (3)  with  approximately  90  percent  of  all  data 
points  falling  within  20  percent  of  the  equation.  A  sensitivity  analysis  shows 
that  this  scatter  is  not  important  for  this  report,  because  flow  rates  are 
weakly  influenced  by  the  value  of  P. 

Note  that  the  laboratory  tests  used  to  develop  equation  (2)  show  some 
ponding  for  high  breakwaters  not  overtopped  and  submerged  breakwaters.  This 
ponding  occurs  because  the  breakwater  reduces  the  wave  height  landward  of  the 
structure  and  the  change  in  wave  height  causes  wave  setup  (Longuet-Higgins , 
1967). 

Overtopping  tests  by  Diskin,  Vajda,  and  Amir  (1970)  show  that  the  net  over¬ 
topping  rate,  qn,  'for  breakwaters  may  be  estimated  from 

*n  -  «io  (i  -  ir)  (4) 


f/hJ, 

Figure  3.  Ponding  levels  for  porous  multilayered  breakwaters  (after 
Diskin,  Vajda,  and  Amir,  1970). 
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2.  Offshore  Breakwaters  with  Gaps. 

Overtopping  of  breakwaters  causes  a  buildup  of  water  in  the  zone  landward 
of  the  structure.  If  there  are  gaps  between  the  breakwaters  and  the  breakwaters 
are  not  connected  to  shore,  the  buildup  of  water  landward  of  the  breakwaters 
causes  flow  through  the  openings  (Fig.  4).  One  method  of  estimating  the  exit 
flow  through  breakwater  gaps  or  inlets  is  to  use  a  combined  continuity-energy 
equation  for  discharge. 


Q  =  VAC  -  Cd  ^2ghb  Ac  (5) 

where  Cd  is  a  discharge  coefficient,  Ac  is  cross-sectional  flow  area  at 
the  water  level  of  interest,  and  V  is  the  mean  velocity  of  water  flowing 
through  the  gaps.  Many  factors  may  influence  the  magnitude  of  Cd  but,  as 
a  first  estimate,  Cd  =  0.8  is  recommended  for  gaps.  The  discharge  coefficient 
for  spaces  between  breakwaters  and  the  shoreline,  Cde,  has  a  recommended 
value  of  Cde  =  1.0. 


Figure  4.  An  offshore  breakwater  system. 
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Note  that  in  this  first  approximation  of  breakwater  gap  flow  that  the  waves 
are  assumed  to  approach  approximately  normal  to  the  breakwaters  and  shoreline, 
so  the  longshore  current  can  be  neglected.  Other  effects  such  as  diffraction, 
refraction,  reflection,  and  wave-current  interactions  have  not  been  considered. 

If  incident  wave  conditions  do  not  vary  rapidly  with  time,  a  condition  will 
be  reached  where  water  flowing  into  the  zone  protected  by  the  breakwaters  will 
equal  the  exit  flow  through  breakwater  gaps  or  inlets.  The  equation  describing 
this  condition  is 


qQ  (x  -  r)  N*  -  b  (cdBds(N  -  1)  +  2CdeAce)  (6) 


where  N  is  the  number  of  breakwaters,  Ace  is  the  cross-sectional  flow  area 
between  the  breakwaters  and  shoreline  at  the  end  of  the  system,  and  B  is  the 
gap  width  between  breakwaters  (Fig.  4),  Solving  equation  (6)  and  putting  into 
dimensionless  form,  the  dimensionless  velocity, 


V 


becomes  a  function  of  a  single  coefficient,  K,  where 


V2FP~[cdBds(N  -  1)  +  2CdcAce; 

qc*N 


(7) 


Figure  5,  which  gives  the  relation  between  dimensionless  velocity  and  K,  shows 
that  any  combination  of  factors  causing  K  to  increase  will  produce  a  smaller 
dimensionless  velocity.  For  example,  keeping  all  other  factors  constant,  if  the 
gap  spacing  B  is  increased,  K.  will  increase  and  V  will  be  reduced. 

K  and  the  resulting  velocity  may  be  easily  solved  using  equation  (7)  and 
Figure  5.  The  computer  program  BWFL0W2  can  also  be  used  to  solve  the  velocity 
and  flow  through  breakwater  gaps  due  to  overtopping.  This  program  is  recom¬ 
mended  if  a  large  number  of  calculations  are  needed.  The  program  is  also  sug¬ 
gested  any  time  irregular  wave  conditions  are  assumed,  because  irregular  wave 
overtopping  rates  are  a  complex  function  of  the  overtopping  rate  given  by  equa¬ 
tion  (1).  The  cost  of  running  BWFL0W2  is  a  few  cents  per  condition  of  interest. 

It  is  recommended  that  V  be  kept  below  0.5  foot  (0.15  meter)  per  second 
for  extreme  design  conditions.  Velocities  much  higher  than  this  value  could 
transport  significant  amounts  of  sediment  out  of  the  breakwater  system  and  may 
cause  scour  around  the  breakwaters.  Recall  that  V  is  an  average  velocity 
through  the  gap  and  that  local  velocities  in  breakwater  vicinity  may  be  con¬ 
siderably  higher. 
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Figure  5.  Dimensionless  velocity  as  a  function  of  K. 


III.  EXAMPLE  PROBLEMS 

Example  problems  are  presented  to  illustrate  the  steps  required  for  manual 
or  computer  computation.  The  examples  indicate  the  relative  importance  of  the 
breakwater  design  variables  on  the  magnitude  of  V. 


***************  EXAMPLE  PROBLEM  1*  ************** 


GIVEN:  Four  rubble-mound  offshore  breakwaters  have  the  design  conditions 
illustrated  in  Figure  6. 

FIND:  V,  assuming  a  monochromatic  wave  height  of  8.0  feet  (2.44  meters)  at 
the  structure. 


SOLUTION :  Manual  computation  is  illustrated  in  this  example. 
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H0‘  =  7.8  ft 
T  =  7  s 

Monochromatic  Wave 


/  =  400  ft 


Ace  =  1,000  ft2 

cde  =  i.o 


M 


B  =  125  ft 


H  =  8  ft 
ds  =  12ft 
h  =  1 3  ft 
tan  0  =  0.667 


Shoreline 


Figure  6.  Design  conditions  for  example  problem  1. 


The  first  step  in  computations  is  to  determine  values  for  the  dimensionless 
parameters: 


and 


7.8 

32.2(7)2 


0.0049 


£s 

Ho 
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7.8 


1.53 


From  Figure  2,  a  and  Q*  are  estimated  as  a  =  0.053  and  Q*  «  0.019.  The 
surf  parameter  and  resulting  runup  are  determined  using  equation  (2) : 


e  = 


y  5 . 


0.667 

_8 _ 

12(72) 


3.74 


and 


_  0.692(8)  (3.74) 

(1  +  0.504 [3.74]) 


7.17  feet 


The  breakwater  freeboard  is  F  =  h  -  ds  =  13  -  12  =  1 .0  foot. 
The  overtopping  rate,  qQ,  given  by  equation  (1)  is 

0.1085 

_ /7 ,17  -  l.\  0.053 

qo  =  /32. 2 (0.019) (7. 8) 3  V7.17  +1./ 


9.53  cubic  feet  per  second  per  foot  of  breakwater  crest. 


From  equation  (3)  the  ponding  level  is 


P  -  7.8(0.6)e-^'^'2  *  l-°/7-8>2  -  3-36  f.et 
From  equation  (7) 

V64.4  3^6  [0.8(125)  12(4  -  1)  ♦  2(1.0)  1.00Q] 

*  “  9.53(400)  4 

or 

K  -  5.4 

The  corresponding  dimensionless  velocity  is 

•  0.18  for  K  -  5.4  from  Figure  5 

and 

V  *  0.18(0.8)  V 64.6  3.36  «  2.11  feet  per  second. 


V 

Cd\2gP 


***************  EXAMPLE  PROBLEM  2 *************** 

GIVEN:  The  same  breakwater  system  as  in  example  1. 

FIND:  V  for  irregular  wave  conditions  with  a  significant  wave  height, 

Hs  -  8.0  feet  (2.44  meters). 

SOLUTION:  The  computer  solution  is  illustrated  in  this  example  with  the 
input  shown  in  Table  1.  Resulting  computer  output  is  given  in  Table  2. 

The  predicted  velocity  is  V  -_1.33  feet  (0.14  meter)  per  second.  This 
velocity  is  lower  than  the  V  obtained  for  monochromatic  waves  in  example  1 
(2.11  feet  per  second),  because  irregular  waves  of  a  given  significant  wave 
height  have  less  overtopping  than  monochromatic  waves. 


Table  1.  Computer  program  input  for  example  problem  2. 


1  4  4  CO.  125.  7.77  8.  7.  12.  K. 


•  8t7  .018  .053  .8  -l.P  If'PC. 


IMIM0MII9SI  sissosoes  OOMMIMMOf  Ml  Ml  MIIIIM  I M  M  M  M  M  M  M  C  M  M I «  M 

I  9  9  •  1 1  t  «  rgR««<;naot«oi  ««»><  vs  uw  ***»<*•  *  it  »  i.*  .<  *  *  w  i*  «  mm 

n  1 1  f  1 1 1 1 1 1  n  1 1  1 1  n  i  ii  1 1  it  1 1 1 1  it  ii » 1 1 1 :  m  itiiiiiii  1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


Table  2.  Computer  program  output  for  example  problem  2. 


••••  iHaEcbl*#  **v£s  •••• 

N  tCftl  «<M1  HO(fT)  H(fT) 
•  *00.00  {25.00  7 . *  ».0 

CD  COE  OCt(ETp)  RU*U*(f  T ) 

.•  1.0  1000. 00  *.1T 


T(itC)  0S<f1>  "*tM)  T**T  0* 
T,0  12.0  l*.0  •<»*7  .01^0 

aotc'oi  wt>  *  v(fPS) 

S.«l  2,2*0  T.12  1.S5 


iLP"* 

,05i0 


***************  EXAMPLE  PROBLEM  3*************** 

G I VEIN :  Three  rubble-mound  offshore  breakwaters  are  located  with  dB  •  12  feet 
(3k  7  meters),  i  ■  250  feet  (76.2  meters),  tan  6  -  0.667,  and  Ace  -  1,440 
square  feet  (134  square  meters). 

FIND:  The  Influence  of  breakwater  freeboard,  Incident  wave  height,  wave  period, 

gap  spacing,  and  shore  attachment  on  V. 

SOLUTION-  The  computer  program  was  used  to  make  these  calculations  and  results 
are  given  it.  Figures  7  to  10.  Figure  7  shows  that  an  increase  in  breakwater 
freeboard  or  gap  to  length  ratio  causes  V  to  decrease.  An  increase  in  in 
cldent  wave  hekght  and  in  wave  period  produces  an  Increase  in  V  as  expected 
(Figs.  8  and  9).  Shore  attachment  with  impermeable  walls  forces  more  flow 
through  breakwater  gaps  than  for  a  detached  system  (Fig.  10).  The  most  dra¬ 
matic  effect  of  shore  attachment  occurs  for  relatively  small  breakwater  gaps. 


*************************************** 

IV.  SUMMARY  AND  CONCLUSIONS 

A  method  is  presented  for  estimating  the  first  approximation  of  the  water 
velocity  and  flow  rate  through  breakwater  gaps  caused  by  wave  overtopping. 
Calculations  can  be  performed  either  by  hand,  using  a  dimensionless  curve,  or 
by  a  computer  program,  BWFL0W2 ,  available  in  the  Corps  of  Engineers  Computer 
Library.  Examples  of  both  calculation  methods  are  given  to  illustrate  the 
relative  influence  of  various  design  parameters  on  the  magnitude  cf  the  gap 
velocity,  V.  It  is  suggested  that  V  not  exceed  0.5  foot  per  second.  High 
values  of  V  may  produce  scour  around  structures  and  transport  sediment  out 
of  the  zone  protected  by  the  breakwaters . 
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